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This paper extends the existing theory of two carrier photorefractivity resonance, which is gener- 
ally apphed to Iron doped Indium Phosphide (InP:Fe), to the case of low non-harmonic illumination. 
The space charge field profile is computed, and the variations of its amplitude, width and position 
are determined as functions of the background intensity. The effect of photorefractive resonance on 
these quantities is evidenced, contributing to the understanding of published experimental results 
in InP:Fe. 

PACS numbers: 42.70.Nq, 42.65.Tg 
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I. INTRODUCTION 

The photorefractive effect in iron doped Indium Phos- 
phide (InP:Fe) has been studied almost two decades ago 
using classical Two and Four Wave Mixing (TWM and 
FWM respectively) experiments [J, These experi- 

ments have shown that InP:Fe exhibits a photorefractive 
resonance at a given intensity, which enhances its TWM 
and FWM gain. This behavior has been successfully 
explained by a band transport two-carrier theory [3, Q 
which is in full agreement with experiments: the reso- 
nance intensity 1^ is the intensity at which the thermal 
electron excitation rate is equal to the optically induced 
hole excitation rate. 

TWM and FWM require a harmonic illumination, 
which allows to describe the induced photorefractive 
space charge field as a complex number: the real part is 
in phase with the illumination while the imaginary part 
is out of phase. In that case, the resonance is character- 
ized by a maximum of imaginary part, corresponding to 
a maximum in the photorefractive TWM gain, and by a 
change of sign of the real part around Ir- 

More recently, photorefractive self-focusing experi- 
ments in InP:Fe have shown both a change of sign of the 
nonlinearity (from self-focusing to self-defocusing, or vice 
versa), along with a shift of the beam[^, 0]. This obser- 
vation has been naturally interpreted through a general- 
ization of the resonance theory, where the local focusing 
effect is assumed to behave as the space charge field real 
part and where the beam shift in interpreted to reveal 
the imaginary part. 

However, even more recent experiments 01, using the 
same setup as in Q but samples prepared in a differ- 
ent way, did not present the inversion of the self-focusing 
behavior or present it at intensity values different from 
the resonance intensity Ir as measured from independent 
TWM experiments 0, Those experiments were suc- 
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FIG. 1: Band transport two-carrier model with one middle 
band trap, characteristic, for instance, of Iron doped Indium 
Phosphide. The arrows shows the electrons motions between 
the energy levels with the associated coefficient as in system 
((!]) and expressions ((2} and (|3}. 



cessfuUy described by a fully numerical model which 
thus does not formally link the inversion to resonance. 

The goal of the present paper is to formally general- 
ize the theory of [3, which we recall in section II, 
to the case of a non-harmonic illumination such as the 
beam which is used in the self-focusing and spatial soli- 
tons experiments[l,[l, Bi- 



ll. THEORY OF PHOTOREFRACTIVE 
RESONANCE FOR TWM IN TWO-CARRIER 
PHOTOREFRACTIVE SYSTEMS 



As gathered from the litterature[10|, a band trans- 
port two-carrier model with a middle band trap can suc- 
cessfully describe photorefractivity in InP:Fe, the middle 
band trap being the Fe^'''/Fe'^''' as shown on figure [TJ If 
reduced to one dimension, this model is described by the 
following equation set: 
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where E is the electric field, x and t the space and 
time variables, e = £^£0 the dielectric permittivity of 
InPiFe, e the elementary electric charge, Nd, Na, n, p 
the densities of donors, acceptors, free electrons, holes, 
respectively, nx, pt and Nt are the densities of iron in 
the form Fe^"*", Fe^+ and total, respectively. j„, jp are the 
current densities, and fJ-p the mobilities, of electrons 
and holes, respectively. K is Boltzmann's constant and T 
the absolute temperature. e„, Cp are the excitation, c„, 
Cp the recombination coefficients of electrons and holes 
respectively. 

The former write as 
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where e^'*, e^J^ are the thermal, and cr„, CTp the optical 
excitation coefficients of electrons and holes respectively. 

The theory of the photorefractive resonance is based on 
the assumption of a low fringe contrast [1, 01 ■ In the case 
of single beam illumination, this hypothesis implies that 
the signal beam lies on a background intensity at least 
an order of magnitude more intense. This is however 
not the case in the previously published experiments!^, 
[1,01 nor in the recent numerical model 0. Furthermore, 
this uniform background beam cannot account for any 
thermal generation, as was done previously for one carrier 
photorefractivity spatial solitons[ll|. 

However, as we will show, this hypothesis is absolutely 
essential if one wants to generalize Picoli's resonance 
theory I, [|. We will thus assume its validity through- 
out this paper. 

The intensity / is written as 



I = Io+he' 



(4) 



Ji being written as Ii = mio, where m is the fringe 
contrast. The low fringe contrast condition thus implies 
the assumption that /i <C /o, or m <C 1. 

The photorefractive space charge electric field is ex- 
panded in the same way, as E = Eq + Eie^^^ + cc + . . 
where E^ is the constant and uniform applied field, and 
the space charge field Eie^^'-^ the linear response to the 
small signal /le''^^. Performing the computation as in 
Refs. U [3 but without any further approximation, we 
retrieve the same expression as in refs [J, 01 for the am- 
plitude of the space charge field Ei , which we recall here: 
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GnOj are given by Eqs. (j2|3p with I ~ Iq, and tito 
and pto are the value of the corresponding densities in 
the dark. 

A first remark is that, as pointed out in the original 
paper i, Eq. ® can be reduced to the same form as in 
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which evidences the resonance intensity 
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only if the mobility fields En end Ep are negligible with 
respect to the diffusion field Ed- 
Notice that the space charge field Ei is proportional to 
the modulation Ji = m/g, while the resonance condition 
involves the background intensity Iq only. The apparent 
relation between the amplitude of the modulation Ii and 
the resonance condition Iq = is an artifact due to the 
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introduction of the fringe contrast m = Ii/Iq. Erroneous 
interpretations may arise if one forgets that the intensity 
of the modulation is measured by m, and not by /q, and 
that m is assumed to be small. Numerical investigation 
of the response to strong modulations, with a fringe con- 
trast m close to 1, can be found in Ref. but nowhere 
else to the best of our knowledge. 



III. THEORY OF PHOTOREFRACTIVE 
RESONANCE FOR A LOCALIZED BEAM 

We consider an illumination pattern Ii (x) which is not 
harmonic any more, and assume that it is small if com- 
pared to a uniform background illumination Iq , as in the 
case of the TWM 0, |1| . The Kuhktarev equations ^ 
can thus be linearized about the solution corresponding 
to the uniform illumination Iq, and solved in the general 
case by means of a Fourier transform, according to 



Eiix) = / Ei{k)e"'''dk 



(9) 



Ei{k) being given by formulas (|5|6|) above, which reduce 
to ([71 [8]) when mobilities are assumed to be large and 
diffusion weak. To let the k dependency appear explicitly, 
let us simplify equations (|6a|6bp by setting Eq — Eg/kri 
and Ed = Eokr2, with 
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Assuming Eq = 10kVcm"\ T = 297 K, uto = 5 ■ 
10^^ cm^^, pto = 6 • 10^^ cm^'^, and = 12.6, which are 



reasonable values for InP:Fe[l, [l^, we get ri = 0.15 fim 
and r2 = 0.026 /im. 

Following, we obtain from ([7]) the following expression 
for Ei{k) and any /q : 
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Approximate expressions can thus be given for a back- 
ground intensity Iq well above, well below and precisely 
at resonance, as the following: 
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Notice that the latter expression is singular for k ~ 0. 

Then, for any given Ii{x), we compute its Fourier 
transform, report it into the above expressions, and com- 
pute the inverse Fourier transform Ei (x) . In general, the 
computation can be performed numerically by means of 
a standard fast Fourier transform (FFT) algorithm. If 
Ii{x) has a Gaussian shape, say for instance 



h = Ae- 



(14) 



the inverse Fourier transform can be evaluated by ap- 
proximately computing the integral by means of the sad- 
dle point method, i.e., for r large enough, by replacing 
the response function Ei/Ii by an expansion in a series 
of powers of k. 
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We get, for /o ^ /,-, up to the second order in k 
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and for Jo ^ Ir in the same conditions: 
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Using the values of ri and r2 mentioned above, the 
first order formulas (i.e. Eqs. (|15|16|1 in which the term 
proportional to 1 /r^ is neglected) are in good agreement 
with the exact numerical solution down for r > 3/im, 
and for r > Ifim if second order formulas are used, see 



Figs. [2]|3l In these figures, as in the following ones, 
the space charge field Ei is normalized by Eu = EoA/Ir. 
The figures are plotted for r = 1/im.For larger values of 
r, the figures have the same shape, but the curves corre- 
sponding to the various computation methods arc close 



FIG. 2: (Color online) The profile of the normalized space 
charge field Ei{x)/Eu, for lo/Ir = 100, and r — Ifim. Solid 
blue : numerical, dashed pink: first saddle point method, 
dash-dotted green:second order one. Dotted black; the inten- 
sity signal — /i (divided by Alo/Ir)- Numerical box width: 
2xm. ~ 40, number of points: 1024. x in fim. 



FIG. 4: (Color online) The profile of the normalized space 
charge field E\{x)/Eu for la ~ Ir, and r — l/im. Solid blue: 
numerical, dash-dotted green: second order saddle point, dot- 
ted black: h/A. Numerical box width: 2xm = 1600, number 
of points: 4096. x in fim. 




FIG. 3: (Color online) Same as fig. [3 but Iq/I,- = 0.01. 
Dotted black line is +I1/A. x in fim. 



together, showing the higher accuracy of the approxima- 
tions. 

Close to the resonance, the singularity in l/k raises a 
difficulty, which can be solved using the fact that 



dEi 
dx 



which yields, through the saddle point method. 
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where erf denotes the error function. 

Compared with the numerical solution (in which 
Ei{k — 0) is set to zero to avoid singularity), the agree- 
ment is reasonable, except that the constant term is lost. 
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FIG. 5: (Color online) The profile of the normalized space 
charge field Ei{x) / Eu- Thin solid dark blue: lo/Ir = 1.5, 
thick solid light blue: 1.1, dash-dotted green: 0.9, dashed 
pink: 0.5. Dotted black line is Ii/A. r — l/im. Numerical, 
with box width: 2xm = 40, number of points: 1024. x in fim. 



and a linear term, which depends on the box size, and 
hence can be considered as a numerical artifact. The 
agreement holds down to r ~ 0.1 /im as can be seen on 
Fig. HI 

For other values of the ratio lo/Ir (i.e. close to 1), no 
approximate analytic expression can be given, but the 
space charge field Ei can be computed numerically using 
equation and a Fast Fourier Transform algorithm. A 
few examples are shown on Fig. [5l evidencing the nature 
of the transition from a positive response of the field to 
a negative one when the background illumination crosses 
the resonance intensity. For r ~ lOfim, the same curve 
shape is obtained for ratios lo/Ir = 1.05, 1.01, 0.99, and 
0.95. 



FIG. 6: The extremal value Eimax/Eu of the normalized 
space charge field vs the normalized background intensity 
lo/Ir- The resonance at /o = Ir clearly appears. 



FIG. 7: The width L (in /im) of the space charge field vs 
the normalized background intensity lo/Ir- The resonance 
at 7o = Ir clearly corresponds to a critical widening of the 
screening. 



/^From the numerical computation of the space charge 
field El , the variations of some characteristics of the lat- 
ter can be computed, and drawn against the background 
intensity Jg . We chose to compute the evolution of those 
characteristics which are relevant for understanding the 
behavior of photorefractive self-focusing and spatial soli- 
tons: maximum space charge field value, its width with 
respect to the beam and its spatial displacement. 

Figure [6] shows the maximum absolute value Eimax of 
El vs Iq. As expected, it is a typical resonance curve, 
showing the photorefractive resonance at /q ^ Ir- 

The effect of the resonance on the spatial extension of 
the space charge field is illustrated on Fig. [71 The width 
L of the space charge field is computed as 



with 
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(20) 



At resonance, /q = /,., the width L diverges, in accor- 
dance with expression (fTS]) and Fig. IH 

The displacement of the space charge field with re- 
spect to the center of the beam is illustrated on Fig. [S] 
The positions Aa;+ of the maximum of Ei , and Aa;_ of 
its minimum are plotted vs Iq. For /q < Ir, the space 
charge field Ei is positive, and A.t+ is the location of its 
maximum, it is positive: the space charge field is shifted 
in the direction of the external field. If on the contrary 
/o < /r, then ii'i < and its location is determined by 
Ax- , which is negative: the space charge field is shifted 
in the opposite direction. For Iq = Ir, the location Aa;_|_ 
of the maximum of Ei is further shifted to the right, but 
it does not represent the maximal value of \Ei \ any more, 
and hence we did not report it on Fig. [Sj For higher r, 
the curves on Figs. [SI [7| and [5] have the same shape, 
except that they go closer to the asymptote. 
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FIG. 8: The positions /S.x+ (upper curve) and Aa;_ (lower 
curve) of the maximum and minimum of the space charge 
field El vs the normalized background intensity lo/E- The 
center of the beam is located at x = 0. {Ax± in /im). 



IV. CONCLUSION 

We have generalized to the general case and particu- 
larly to the case of a localized beam the theory of the pho- 
torefractive resonance in two-carrier photorefractive sys- 
tems such as Iron doped Indium Phosphide [1, 01, which 
pertains only to two-wave and four-wave mixing experi- 
ments and the harmonic illumination they produce. 

As could be inferred from the the harmonic theoryj^, 
3, the photorefractive space charge field indeed changes 
its sign around resonance. The amplitude, width 
and location of the space charge field show a singularity 
and/or an inversion of their behavior around the reso- 
nance. However, and in contrast to previously thought, 
the resonance condition docs not compare the resonance 
intensity to the intensity of the localized signal itself, but 
to the intensity of a uniform background illumination, 
which is itself not so far from the mean intensity of the 
overall illumination of the sample. 

Furthermore, this background intensity is assumed to 
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be large with respect to the signal one. This shows that, 
in the experiments in which only a localized illumination 
is present {i.e. a single Gaussian beam), Picoli's theory of 
the photorefractive resonance does not apply. Hence 
there is a mere analogy between the phenomenon of in- 
version of the photorefractive response observed in Refs. 

Q and the photorefractive resonance stricto sensu, as 
no background illumination is present. Especially, there 
is no particular reason why the observed inversion in- 
tensity should coincide with the resonance intensity 
measured by TWM and computed from the resonance 
theory. 



and hence g and h should be proportional, while they 
are arbitrary independent functions. Consequently, the 
assumption p{x) oc n{x) in Ref. is by no means 

justified. 

Further, even if the property p{x) oc n(x) is admitted, 
a subsequent mistake is found in the reasoning. Assume 
indeed that this statement is satisfied. Combining Eqs. 
(jldlllfl) . looking for a stationary solution and neglecting 
diffusion, we obtain 



^ Un+jp) = 0, 



(A4) 



APPENDIX A: DISCUSSION OF A PREVIOUS 
THEORY 



and consequently 



Reference [ij] pretends to give a theory of the photore- 
fractive resonance in the case of a localized beam. How- 
ever, this theory is erroneous, and we intend to prove it 
in this appendix. 

A mathematical error is found in the Appendix A of 
Ref. . The authors indeed state that 'These equations 
share the same solution (...); therefore they must have the 
same coefficients', this argument is false. 

Let us write a few elementary mathematics to justify 
this formally. Denote by /, g and h three arbitrary func- 
tions of X. Then set ip = d{fg)/dx, ~ d {fh)/dx. 
Straightforward computation shows that 



d£_ _^ dg/dx j _ ^ 
dx g g ' 



(Al) 



— {nE) = 0, (A5) 

as in Ref. [T^ . Let us denote by , and E^o the val- 
ues of n, p and E at infinity in x, we get n = n^Eoc/ E, 
and p = paoEao/ E. Reporting these values into system 
([1]); and taking into account the fact that jn and jp do 
not depend on x because their sum does not (Eq. (jA4j ) 
and they are proportional due the assumption, Eq. (jldp 
and pe[) yield respectively the two equations 
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(A6) 



and 



df ^ dh/dx ^ Tp 
dx h ' h 



(A2) 



The two equations (|A1IA2[) share the solution /, hence if 
the 'property' involved in the cited paper were true, we 
should have 



dg/dx dh/dx 



(A3) 
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which arc not compatible, since E only depends on x. 

In short: system ([T]) does not admit any non- uniform 
solution with p(x) oc n(x). Consequently, Eq. (2) of Ref. 
[T^ is incorrect, and the consequences drawn from it are 
not founded. 
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